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PREFACE 


The work described in this report was performed by the Space 
Sciences Division of the Jet Propulsion Laboratory. 
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ABSTRACT 


In this report, we present the method used to calculate the velocity, 
temperature, and density of the solar wind plasma from spectra obtained 
by attitude- stabilized plasma detectors on the earth satellite OGO 5. The 
method, which uses expansions in terms of Hermite polynomials, is very 
inexpensive to implement on an electronic computer compared to the least- 
squares and other iterative methods often used for similar problems in 
the past. 
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COMPUTATION OF SOLAR WIND PARAMETERS 
FROM THE OGO-5 PLASMA SPECTROMETER DATA 
USING HERMITE POLYNOMIALS 

DESCRIPTION OF THE INSTRUMENTS 

Ludwig (1963) has given a description of the OGO satellites. The 
plasma experiment reported here consists of two nearly identical packages: 
one looking toward the sun, and the other looking radially away from the 
earth. Each package contains both a curved-plate electrostatic analyzer 
and an ac Faraday cup . 

Let us consider only the solar-oriented package. The basic concepts 
of its operation can be understood with the aid of the flow diagram in 
Fig. 1. 

The Faraday cup makes high time resolution measurements of the total 
charge flux of those positive ions whose velocity component parallel to the 
satellite-sun line corresponds to an energy per unit charge (E/Q) between 
100 and 11,000 volts; this range of E/Q usually includes both the protons 
and the a particles in the solar wind. 

The average direction of positive-ion flow can be computed from the 
relative currents reaching each of four collectors in the Faraday cup. The 
details of this calculation are given in Appendix 1. 
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At the same time that the ions' total charge flux and direction are 
measured by the Faraday cup, a series of voltages is applied to the elec- 
trodes of the electrostatic analyzer to obtain the E/Q spectrum either of 
positive ions or of electrons. This curved plate analyzer has much in 
common with one flown on Mariner 2 (Neugebauer and Snyder, 1966; Graham 
and Vescelus, 1967). Positive ions entering the detector pass through 
(1) an entrance aperture at the spacecraft ground potential, (2) the 
cylindrical deflection electrodes, which are 120° in length and which have 
the ratio of the radius of curvature of the outer electrode to that of the 
inner electrode equal to 1.134, (3) an exit aperture at ground potential, 
and (4) a suppressor grid at a potential of — 12 volts with a transparency 
of 927,. The ions then enter a collector cup, and the current to this 
collector is determined by the electrometer system for each voltage setting 
(channel) on the deflection electrodes. 

There are 128 overlapping E/Q channels equally spaced on a 
logarithmic scale between 2.54 volts and 16.9 kv; a particles should 
appear 10 channels above protons with the same velocity. A few more 
details about the operation of the curved plate analyzer are given in 
Appendix 2 . 

The assumption that each ion species in the plasma can be described 
by a Boltzmann distribution (temperature T) in a reference frame moving 
relative to the detector with velocity v q in a direction given by 9 and 
a Hows one to calculate values of T, v , and density n, from a spec- 
trum. It is then possible to apply a temperature correction to the 
calculation of 9 q and cp^ from the Faraday cup currents and then repeat 


JPL Technical Memorandum 33-519 



the calculation of v , T, and n (see Fig. 1); however, for the ratio of 
thermal velocity to flow velocity usually found in the solar wind, this 
correction is negligible compared with other uncertainties in these angles. 

Next, the apparent, or measured, total charge flux can be corrected 
for angle of incidence and temperature to obtain an estimate of the true 
total charge flux = ^<l i n i v i > where the sum is taken over all species of 
positive ions. 

It is known that all species of ions in the solar wind have approxi- 
mately the same velocity (Bame et al., 1968; Ogilvie et al., 1968) and that 
the solar wind is electrically neutral (Montgomery et al., 1968). Thus, 
(£q^n^v\)/v = n^, the electron number density. 

The electron number density can also be calculated by summing the 
positive-ion q^n^ as determined from the curved-plate analyzer spectra. 
However, a large thermal anisotropy, which has often been observed in the 
solar wind (Scarf et al., 1967; Hundhausen et al., 1967) will introduce 
a greater error into this method than into the method using the total 
charge-flux measurements. 


METHOD OF COMPUTING PLASMA. PARAMETERS FROM CURVED PLATE ANALYZER SPECTRA 
The current <3 measured by the electrometer of the curved plate 
analyzer when the voltage across the plates is set for channel I is given 
by 


«*(!) 


qA ///- (9) 8 =(^' cp ) 


v cos 9 cos cp 


( 1 ) 


f(v,e,cp) d 3 v 
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where : 


q = charge of the incoming ion 
A = effective detector aperture area = 2.26 cm 2 
v = velocity of the ion 
v = velocity at center of channel I 

The angles 9 and cp are defined as shown in Fig. 2. In this 
coordinate system 

d 3 v = v 2 cos 9 d9 dcp dv (2) 


The product g x (9) g 2 (v/v^.,cp) is the response function of the analyzer. 

gl (9) = 1 + 9/10.1° for 9^0° 

= 1 - 9/10.1° for 9 ^ 0° 

l I o (3) 

= 0 for |0| > 10.1° 

The function g^v/v^cp) cos cp was determined experimentally by calibration 
with a cold argon ion beam of large cross-section. This function is 
plotted in Fig. 3. 

f(v,9,cp) is the distribution function of the ions. If we assume 
the particles have an isotropic Boltzmann distribution superimposed on a 
convection velocity of magnitude v q and direction given by 9 q and cp , 

3 


f(v,e,cp) 



v-v | 2 /2kT 
o 1 


(4) 


v-v 2 = v 2 + v 2 - 2v v [cos 9 cos cp cos 9 cos cp 


+ cos 9 sin cp cos 9 sin cp + sin 9 sin 9 ] (5) 

o o o 
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When Eqs . (2) through (5) are substituted into Eq. (1), the 9 
integral can be done either analytically (if it is assumed that 9 is a 
small angle) or numerically (for any 9). The 9 and v integrals have been 
done numerically. The result is 

c9(l) = qA n v q G[v(I) ,9 o> 9 o ,S] 


where 


v(I) = log (v I /v o ) 


5 = 



v = velocity at center of channel I 

v = convection velocity 
o 


Examination of many computed values of G suggested it could be rewritten 


as 

G[v(I),9 o ,9 0 ,5] = }t[v(I), 9 o ,|] = G[v(I),0,9 o ,|] j X 10 _9 o Gl(5) 


A least-squares fit of the function 


1 



o 


G[v(I) ,9 o ,9 0> |] 
108 G[v(I),0,9 o ,SJ 


to a power series in § gives 

G x (§) = (5.05006 - 10.6985 % + 8.62214 f - 1.506114 f) X 10' 3 
for 9 q in degrees. 

T[v(I),9 o ,?] has been computed and tabulated as TAB(I,J,L): 

I = 1 to 16 for v(I) = [- 8 log = - . 1204l] to [+ 7 log 2’ G ^ = 

.10536] in steps of Av(I) = log 2’° 5 = .015052 = 1 channel. 
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J - 1 to 11 for cp — - 10° to + 10° in steps of Acp = 2° 
o o 

L = 1 to 10 for ? = 1.2 to 3.0 in steps of A? = 0.2 


Appendix 3 is a listing of TAB(I,J,L). [Note: A minor goof 

resulted in each value of the array TAB(I,J,L) being too big by a factor 

of 2. This was taken into account when calculating density. J 

In computing the plasma parameters corresponding to a measured 

series of =9(1) (i.e., a spectrum), it is assumed that the angles 0 and 

o 

known from the Faraday cup measurements. A computer program was 
written to find the values of n, v q , and § = log( mv^/2kT) which minimized 
the function 



<^(I) - qA n v q 


10 





3 


The computations were iterative and time consuming. Thus we wrote a 
second computer program which computes n, v o> and g by comparing the 
coefficients of the Hermite polynomials of the measured and the theore- 
tical spectra; it is a relatively quick, non-iterative calculation. 


HERMITE POLYNOMIALS 

A function w(x) which is shaped sort of like a Gaussian, can be 
expanded in a series of Hermite polynomials as: 


w(x) 


C o e -(x-x) s /2o 2 

a yjzii 




where the Hermite polynomials of order 3 and 4 are 
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H 3 (x) = x 3 - 3x 
H 4 (x) = x 4 - 6x 2 + 3 

and in general H n+1 (x) = x H^Cx) - H^Cx) 

= x H n (x) - n H n _ x (x) 

The low order coefficients can be expressed in the following form: 


C 

o 



dx 


C x = x 


C 

o 



x w(x) 


dx 


C 3 = a 2 


J x 2 w(x) dx 



o 


C, = 


1 


O 


w(x) dx - 3^02 - Cj 


c 4 = — 


- C 


x 4 w(x) dx - 40x01 C 3 - 6cfc 2 - 


- 3 


DISCUSSION OF THE METHOD 


The method used to compute the plasma parameters via Hermite poly- 
nomials consists of: 

1. Using the measured spectrum to compute the Hermite coefficients 
C™ through C 4 m , where the superscript m refers to measured values. 
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2. For the angle of incidence determined by the Faraday cup, 
computing a set of coefficients through for each of the 10 values 
of the parameter 5 in the comparison table, where the superscript t 
refers to theoretical values. 

3. Finding the value of 5 for which Cg = Cg . 

4. Using § from step 3, C™, C™, C^(S), and C^Cl) to directly 
compute v , T, and n. 

5. Using C o through C™ to fit a curve through the measured spec- 
trum and then compute the relative variance of the measurements from the 
curve . 

6. Computing the differences C™ - CgCl) and c“ - 0^(5) as measures 
of the excess skew and excess kurtosis, respectively, to tell in what 
directions the measured spectra differ from spectra corresponding to con- 
vected Boltzmann distributions. 

The principal complication, which is glossed over in the above out- 
line, arises because the Hermite polynomials are orthogonal over the 
interval -°° to +=°, whereas the measured data cover a limited range of 
energy which often does not even include all important contributions to 
the peak. On 0G0 5, the number of points /spectrum were limited by 
(1) only allowing 4 channels above the peak channel, to avoid contamina- 
tion by heavier ions, (2) the finite noise level of the electrometer, and 
(3) at high bit rates, the first current measurements following a switch 
from the low to high sensitivity scales of the electrometer could not be 
used. We get the best results when the energy range used to compute the 
Cq through C 4 in step 2 is truncated to match the energy range of the 
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measured spectrum as closely as possible, rather than using the entire 
energy range over which TAB(I,J,L) > 0.0000. 


COMPUTATION OF COEFFICIENTS OF MEASURED SPECTRUM 

A spectrum consists of current measurements in 1 to 8 channels 
below the channel with the largest current (the "peak channel") plus 1 
to 4 channels above the peak channel; no more channels are taken above 
the peak channel to avoid contamination by heavier ions (for OGO-5, the 
alpha peak is expected to be 10 channels above the proton peak) . To 
compute the Hermite coefficients for a measured spectrum, the independent 
parameter x is chosen to be the difference in channel number: x(I) = 

(channel number of measurement I) - (peak channel number), and the func- 
tion w is the measured current ^(1) . Then the first five Hermite 
coefficients are calculated as follows: 


C™ = J' w(x) dx ss ^ ^ c9( I) * A(I) 


where the A(I) = 1, 2, or 4, according to Simpson's rule of numerical 
integration. 


C i =J t Y] 0(1) * X(I) * A(I) ] = 

r i 


x 


o I 


= ^ X ) [c9(i) * x(i)S * a(i)] ■ (c ?> 2 = ct 


C? = 


(c“) § 


~]CO(I) * X(I) 3 * A(I)] - 3c“ C™ - (c”) 3 l 

_ C o 
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cl = 


(C”) 2 


i<hd * x(D 4 * a(i) ] - 4^ 


* (C™) 1 * C” 


o I 


- 6(C“) 2 * C™ - (C?) 4 


COMPUTATION OF COEFFICIENTS FOR COMPARISON WITH MEASURED VALUES 


First, note that the pos ition and shape of the spectra depend only 
on the function T(v.j./v o ,cp o ,§) , represented in the computer by TAB(I,J,L). 

_ q2 q / p \ 

The factor qAnv^ 10 o 1 , on which the amplitude of the spectra depend, 

will be worried about later. Thus, the first step in finding the theore- 
tical coefficients to compare with the C m is to interpolate TAB(I,J,L) to 
find a two-dimensional table TA(I,L) for the value of cp Q determined by the 
Faraday cup. Then, for each value of the unknown parameter | (index L in 
the table), we compute the C^(L) as follows: 

1. Using the measured currents in the peak channel and the channels 
on either side of the peak, interpolate to find the channel number of the 
peak of the measured spectrum. If c9(max) is the largest current measurement 


x(peak) m 


^(max+ 1) - c9(max - 1) 

2[2cS(max) - c5(max + 1) - <5( max - 1)1 


2. Similarly, find the peaks of the theoretical spectra in the 
comparison table: 


x(peak ,L) 11 = -= TA(,LP+1 ? L) - TA(LP-1,L) 

2[2TA(LP,L) - TA(LP+1,L) - TA(LP-1,L)J 

where LP is the value of I for which TA(I,L) is greatest. 

3. Then, the values of x ± which we wish to use in computing the 
Hermite coefficients are given by 
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- LP + xCpeakjL) 1 ” - max - x(peak) m + i 

4. Interpolate the TA(I,L) to find its values at each x : call 

1 

these values TH(IX,L) . 

5. C (L) = , TH(IX.L) * A(IX) where the A(IX) are again the 

IX 

Simpson's rule coefficients. 

The equations for the higher order theoretical coefficients C^(L) 
through C 4 (L) are exactly like those for the measured coefficients c m , 
except that TH(IX,L) is used in place of c9(I) , with the summation being 
over the IX. 

CALCULATION OF VELOCITY AND TEMPERATURE 

Next, we perform an inverse interpolation to find the value of | 

t ni t 

at which C 2 (L) equals C 2 . Then, interpolate the C-^L) to find the value 

of Cj corresponding to this §. The relation between C™, C^d), and the 

channel number corresponding to the flow velocity v q is perhaps made 

clearer by the diagram in Fig. 4, together with reminders that the TA(I,L) 

were defined such that v /v q = 1 when 1=9, that the maximum values of 

TA(I,L) are at I = LP (LP usually = 8, 9, or 10, except when the angle is 

very large), and that C x is the x or center of the measured spectrum rela- 

• Vc 

tive to the observed peak channel. From Fig. 4 we see that x = (channel 
number corresponding to v ) = [peak channel number + cf - C^d) - [LP + 
x(peak,L) ]| + x(peak) - 9} where the subscript ? implies interpolation 
of the function within [ ] to the value 5. 

For 0G0 5, the conversion of x to velocity in km/sec is given by: 
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v = 10 1. 32804 + 0.015052(x*-A) 
o 

where A = 0 for protons and 10 for alpha particles. Temperature in °K is 
given by 

T = 60.5 v 2 M/10^ 
o 

where M = 1 for protons and 4 for alphas. 

CALCULATION OF DENSITY 

To find density, interpolate the C^(L) to find the value correspond- 
ing to ?. Then, 

1.08 * C m * 10" 9 o 
o 

n = 

q * c z o a) * v o 

2K 5 

where Q = 1 for protons and 2 for alphas. The constant 1.08 = — - X 10 
where the factor K = 10 -14 /.512 converts the measured currents from 
telemetry digits into amperes (all C m were computed from the digital 
output of the telemetry system, which is directly proportional to the 
collected currents); e = 1.6 X 10 -19 amp; A is the effective detector 
area = 2.26 cm 2 ; the factor lcf converts v q in km/sec to cm/sec; and the 
factor 2 corrects for the error in TAB(I,J,L) mentioned earlier. 

GOODNESS OF FIT PARAMETERS 

The following three parameters are computed to measure the good- 
ness of fit of the convected Boltzmann model to the measured data: 

1. The relative variance of the data from the polynomial fit 
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where I = channel number and I = peak channel number. 

P 

2. The excess skew = C^ 1 - G^C?) . 

3. The excess kurtosis = C™ - 

Excess skew can arise when the actual distribution function rises 
more rapidly and decays more slowly with increasing energy than a convected 
Boltzmann distribution; i.e., a high energy tail would give a positive 
value of excess skew. Figure 5 shows the effect of a large amount of 
excess skew (ESKEW) for a typical (but hypothetical) spectrum with cp Q = 

4°, 5 = 2.0, and with three points each above and below the peak channel. 

The spectral distortions due to excess kurtosis (EKURT) , or curva- 
ture, are shown in Fig. 6. 

ACCURACY OF THE METHOD 

As a check on the accuracy of the method, we have fit a selection 
of 28 spectra by both methods — least squares (LS) and Hermite polynomials 
(HP). The results of these computations are given in Table 1. The values 
of v, T, and n given are those determined by the least-squares fit. The 
Av, AT, and An are the Hermite polynomial values minus the least square 
values. The LS error is defined somewhat differently than is the HP error, 
but to a first approximation they both represent the variance of the data 
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from the curves. The HP error is always exaggerated because the truncated 
energy range used in computing the Hermite coefficients leads to a poor 
fit (underestimation) at the fringes of the spectrum: this effect does 

not necessarily lead to poor estimates of the plasma parameters because 
of the pains taken to truncate the theoretical spectra to match the measured 
spectra. This effect is illustrated by the fit of spectrum #1, shown in 
Fig. 7. For this spectrum, Av/v = 0, AT/T = 0.007, and An/n = 0.024; i.e„, 
there was very close agreement between the LS and the HP determined param- 
eters despite the apparently poor fit of the HP curve. 

In comparing the fits of the other spectra, we note that |Av/v| was 
less than 0.02 for all cases and was less than 0.01 for all but one spectrum. 
The agreement of the temperatures computed by the two methods is poorer 
than for the velocity. Of the 28 spectra, five had | AT/T | > 0.20. One of 
these, spectrum #22, for which AT/T = 0.230, was only 3 channels wide and 
had a huge HP error. The HP method is not well suited to 3-point spectra. 
The other four spectra are plotted in Figs. 8a through 8d; they were all 
weird shaped spectra. In the case of spectrum #18, the HP values of 
velocity and temperature are probably more realistic than the LS values. 

A large value of | AT/T | can arise from a poor fit with either the 
LS method, the HP method, or both. In Figs. 9 and 10, | AT /T | is plotted 
vs. the LS error and the HP error, respectively. Both plots show some 
correlation between | AT/T | and error, which confirms again that for smooth, 
well shaped spectra, the HP method gives satisfactory results. 

The average values of |Av/v|, | AT/T | , and |An/n| for the 28 spectra 
are 0.003, 0.109, and 0.062, respectively. However, the 28 spectra were 
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selected so as to get a wide variety of different spectra; a random 
sampling would probably give appreciably smaller average differences. 

COMPUTER PROGRAM 

The Fortran subroutines used to perform the HP fit are given in 
Appendix 4, together with commentary translating the program to the nota- 
tion of the text. 
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APPENDIX 1 


FARADAY CUP 

A scale drawing of the Faraday cup detector is shown on the left 
side of Fig. 11. The area of the circular aperture is 20.3 cm 3 . The 
structures labeled G, S, and M are ground, suppressor, and modulator grids, 
respectively. The thick ground grid is actually a solid plate with holes 
in it. The over-all transparency of the combination of grids was experi- 
mentally determined to be 217o for normal incidence, and the relative 
effective aperture decreases almost linearly with the angle of incidence 
such that the half-width at half-maximum is 20.3° . 

The structure denoted by C lj2j3 is actually three separate collectors, 
each in the shape of a 120° sector. Cj is a separate collector for parti- 
cles that arrive at angles greater than ~ 12° . The right side of Fig. 11 
shows these four collectors as viewed from above, together with a smaller 
off-center circle that represents an image of the entrance aperture for 
particles incident from a direction defined by the angles a and P (R is a 
constant) . The sizes of the shaded areas Aj , A 2 , A 3 , and Ai are propor- 
tional to the currents to collectors 1 through 4, respectively. For 
observations in the solar wind, A 4 is usually negligible compared with A x , 

Ag , and A 3 . 

Each collector individually and all collectors combined are sequen- 
tially connected to an electrometer system whose output is telemetered as 
a 9-bit number proportional to the current detected, along with one more 
bit to indicate the sensitivity scale of the electrometer. Measurements 
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and 


l —1 2 

of 1 telemetry unit are equivalent to currents of (0.512) - X 10" 

(0.512) -1 X lO -10 amp on the high- and low-sensitivity scales, respectively; 
the instrument switches to the proper scale automatically. 

The responses of the Faraday cup and the curved-plate analyzer were 
experimentally checked with ion, electron, and neutralized ion beams. The 
following set of equations (relating the ratios L = A. 3 /A-i and M = A 2 /A-i to 
the angles a and P) was derived from the calibration measurements, using 
both 1- and 5-keV argon ion beams with energy spreads of ~ 3 eV full width 
at ha If -maximum: 

P(deg) = (10.2 - 10.2 L) if L s 0.5 


= (14.7 - 19.2 L) if 0.25 ^ L ^ 0.5 


= (18.9 - 36.0 L) if L < 0.25 


«<“**> * (iVlXm+l) x 60 


Figure 12 illustrates the accuracy to which o' and (3 can be deter- 
mined for a particular example in which A x + A 2 + A 3 = 40 telemetry units, 
which is equivalent to a fairly weak solar wind with a total charge flux 
of 1.35 X 10 s charges cm“ 3 sec _1 . The distance between points in this figure 
indicates the uncertainty in a and P arising from the digitization of the 
collector currents. 

The angle a is measured counterclockwise (as viewed from the sun) 
from the center of collector 1. The orientation of the Faraday cup is 
such that the angles 9 q and cp o used in analysis of the curved plate 
analyzer data are given by: 
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| 9 q | = 3 | sin(ar +45° ) | 


cp =3 cos (o' + 45° ) 
o 


The sign of 9 o is not important due to the symmetry of the curved plate 
analyzer . 


0 
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APPENDIX 2 


FURTHER DETAILS OF THE CURVED PLATE ANALYZER 

The electrometer connected to the curved plate analyzer has three 
sensitivity scales with automatic switching between scales. The input 
current to the electrometer is directly proportional to the output volt- 
age, which is digitized as a 9-bit number, along with 1 bit for polarity 
and 2 bits for scale indication. Table 2 gives the sensitivity, range, 
and rise time of each of the three electrometer scales. The data pre- 
sented in this report are on the middle and the most sensitive scales. 
Owing to the electronic transients generated when a scale change is made, 
the first current measurement made after a switch from the most sensitive 
to the middle scale and the current measurements in the first m channels 
made after a switch from the middle to the most sensitive scale have been 
eliminated before analysis of the spectra, where m = 1, 4 and 8 for space- 
craft data rates of 1, 8, and 64 kbit/sec, respectively. At 64 kbit/sec, 
when the peak of the proton spectrum is on the middle scale, which it 
often is, the elimination of eight channels' worth of data upon switching 
back to the most sensitive scale usually results in the elimination of 
most of the cy-particle spectral peak. 

A correction must also be made to the zero-current level of the 
measured spectra to account for an anomalous "photodip" caused by photo- 
electrons entering the electrometer system through an exposed slit in the 
electronics' cover plate. This photodip is about 20 channels wide and 
reaches a peak negative current of between 2 and 8 X 10 -12 amp at or near 
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channel 72 (E/Q - 348 volts) . It varies somewhat in both amplitude and 
channel. A few examples are shown in Fig. 13. Fortunately, the photodip 


can usually be separated from the posit; 


peaks of the solar wind ions ; 


the separation process is aided by observations made when little or no 


positive-ion flux enters the analyzer, as often occurs in the magnetosheath 
and magnetosphere . 


The positive-ion data reported here consist of "narrow-fine" spectra 
for which the instrument measured the currents in only the 32 consecutive 
channels surrounding the channel with the largest current, as determined by 
the previous spectrum. It is also possible to sample all 128 channels con- 
secutively ("wide-fine" spectra) or to sample only every fourth channel 
( wide-coarse spectra) . Compared with the wide-fine mode, the peak-seeking 
feature of the narrow- fine mode of operation allows a fourfold increase in 
the rate of obtaining spectra, with no loss of energy resolution. Unfor- 
tunately, the instrument will choose the anomalous photodip as the "peak" 
if it has greater magnitude (sign is not examined) than the true ion peak. 
Thus, narrow-fine spectra occasionally miss the solar wind peak entirely. 

Furthermore, wide-fine, wide-coarse, or narrow-fine electron spectra 
can be obtained by reversing the polarity of the deflection voltage and 
grounding the suppressor grid when low-energy electrons are detected. The 
methods described in this paper could be applied to electron data, but 
other methods (which will be described in a separate report) are also avail- 
able for cases when it can be assumed that the electron distribution function 
is isotropic in the satellite frame of reference. 
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.00000 
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.00000 
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.0 31 10 

.08074 

. 1 5 7 5 4 

.22 90 5 

1 

2.2 

6 . 

. 239 GG 

.17861 

.09221 

.03217 

.007 36 

.00107 
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.00021 

.00117 
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.00015 
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2 
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- 8 . 
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1 
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2 
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.00000 
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1 

2.0 
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.00024 
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2.6 

8 . 

.00000 
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.00070 

.0091 1 

.06219 

.21622 

.3 72 91 

1 

2.6 

6 . 

. 71 21 1 

. 1 250 3 

.02282 

.00180 

.00000 
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.00000 
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•5 
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.00000 
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.00085 

.00885 

.02238 

.05708 

.05685 

1 

2.6 
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.00553 
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.00002 

.00000 

.00000 

.00000 

.00000 

n 

2.9 

- TO . 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00010 

.00183 

! 

2.9 

- 10 . 

.01132 

.02165 

.01085 

.001 17 

.00002 
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.00000 

.00000 
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.0 

.00008 

.00139 

. 01 917 

1 
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.00009 
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.00000 
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.00000 
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2.8 
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.0001 9 

.00000 
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.00000 
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.00000 

.00000 

.00000 
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I 
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APPENDIX 4. FORTRAN SUBROUTINES 


3F OR .IS FIT. FIT 

SUBROUTINE c TT(U t ’ MAG, VEL. TEMP. DEN, NUTS) 

COMMON /MOMENT/ ICO . ISTOP , TA I 1 G . 1C ) . J »MA X T I 10 ) 

COMMON /MAXCH/MAXCH 

COMMON /TABLE/ TAB < 1 G , 1 1 . 10 ) 

COMMON /IPART/ 1° AP T 
1 /AnGC/ THE t A , PHI, TFLX . CONE .CLOCK 

1 /NOP/ NCHAN 

1 / 1 C H A N / T C H A N ( ? 0 1 

DIMENSION T 3 ( 10! , T N I 10) 

DIMENSION I/A (IS ) , TMUGl.XSAtfE I 10) 


a fwa ex')-- 

■>= MovacaJ’ yv-tAAi. 


= 2 - cufpKa s 
rheTA = 0 o 
PHr-f 0 


H C-HA-I/ " "VUvnUU*. X'A- 

TC.H Ai 4 (jO ' 


DIMENSION UPMAOt 1 I 

DIMENSION V A R I 1 7 ) . TZI10 ) .TF 1101. TStlOl 

J M - IF I X ( PH T ♦ 1 2. ) /2- 1 


TF(JM.L T .1> JM=I 
IFtJM.ST.8) JM c 8 


00 32 1= 1 . 10 
T M A X -C « 

DO 31 I r 1 . 1 G 
00 3 C 0 = 1.11 

30 l/» R C J J = T A9 ( I , J.Ll 

CALL INTRP( (PHI+ 12, ) / 2 • . tf AP , T A ! I.L ) « JM, 8) 
IF ( TA ( T, L ) .L T . TMAX ) GO TO 3! 

IMA X:T A ( I ,L ) 

M® X T ( L ! = I 

31 CONTINUF 

32 CONTINUE 
EMA X = P . 


T NT (?-P D> a 
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DO 3M r-J.NCMJN 


IFIUPK*G CII .LE .EMAXI GO TO 3M 
EMA X=UPMAG ( I 1 
K H X E - I 

MAXCHrICHANtT) 

3<l CONTINUE 

E XPK = IUPMAGt MAXE* ll-UPM AGtM AXE-11 l 
I (MAXE-1 11 
DO 33 L -1 .10 
LP-MAXT(L) 

TXPK=(TACLP+1iL1-T#«lP-1«L11/2./12 
1 ) 

00 36 I -1 . 1 6 

36 I/A (Il-TQ (TiLI 
DO 37 Mrl.NCMAN 
X:TXFK-EXPK-. c LOAT(LP.M-MAXEi 
IF (M.EQ.MAXE1 XSAVE(LI=X 
IG0SX-1. 

IFtIGO.LT.il 100=1 
IFtIG0.GT.131 TOO" 1 3 
CALL INTRPt X. VA. TH(M) .IPO.A) 

IFtTHtM) .LT.O.) TH(M)=0. 

37 CONTINUF 
TOIL > = SKP ( TH . 1 « NC 4 A N ) 

T c t L I =F MOF t 1 . TH 1 / T Zt L ) 
TE(Ll=FMOE(?»TH)/TZtL)-TFtLI»»Z 
IF (TStL 1.LE.0. 1 GO TO 4 3 
T3tL)=tT5(Ll»*t-1.5)l»tFH0Et3.THl/TZILI 


E HAX = <XKl) 

uivc , T £ Cr) - max, 

flHAU ' 0 U 

MAXC.V\ * siuA/nJ*** ^ d)(.T) 

./ I Z •» EMAX-UPMAGt KAXE ♦ 1 l-UPMAG 

EXPk - XL peoAt) 1 * 1 

TAlLP.L)-TAtLP*l.l)-TAILP-l.LI 

T*PK = Ol( pe*J? L)^ 


''-'ll 


T* U) ^ C/(L) 

TF C l) = C^Cl) 

TS Cl) : C/U) . 

T3Cl) - C 3 t U) 

- 3.»TF 111 • TS «L1- TFtLl *»3 I 
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T*ltL» = fTStL I -2 * »*tFMOEI«l*THl/TZil.* -*l . *TF IL»*I TSIL (1.51 I *T3IL 
1 1 -5. ♦ TF f'_ I * ♦ 2 * TOf L 1- T F IL 1 * *9 >- 3. C ^ L ) 

WRITE fE* 90 1 L.T7fL).TFtLI.TSlLliT?fLl.T9IL> 

9C FORM Uteri . 5E1 0. 3 .141 
39 CONTINUE 

F , =jMPIUPMAG*1#NCH#N) 

EFrFMOE ( ! .UPMAG 1 /EZ 
F S=FMO c f Z. UPMAG) /E Z-EF* FF 
IFfES.GT.O. )G0 TO 99 
WRITE ( 5 1 1 9 9 ) 

199 F0RMATfl3H NEGATIVE EE.) 

NUTS:" 

RETURN 


ET-55 ’ 

o 


c; 


n od 



C ” 


r\ju.T5 - I, bJd -w O.K- 


P tunJ&yi 't 


) 

99 E3=SS**f-1.6l*<FMOEt 3 i UPM A 3 1 /SZ-3.«EF*E5-EF**3) 

E9:ES**t-")*IFMCE(9.UPMA3I/EZ-9.*EF*l£S**I.5)*E3-6.*IEF*»?)*F5-EF* 

1*91-3. 

J~ 0 

URITEfS .901 J.EZ. EF ,ES .S3.E9 
KNT:0 

DO 52 J : 1 . 9 

IFfTSC J) .GE .ES.AND.TSI J+U.LT .ESI GO TO 51 


GO TO 52 

5 1 J?=J 
KNT:KNT+1 

52 CONTINUE 

IF (KNT. FQ. 1) so TO GO 
IFfKNT.E3.0l WRXTF(G.53) 

3 FORMAT(35H ES NOT WITHIN RANGE OF TS TABLE. | 
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IFIKNT.GT.il WRT TT( R , j 


M FORMAT ( 2 3 H T 5 TABLE MUL r I - Y AL UE D . ) 

55 NUTS -2 
RE TURN 

49 WRITE (G .50 J 

50 FORMA T I ? 2H NEGATIVE T5tJ> VALUE.) 

GO TO 55 

50 IG0=J7 

IFIIG0.GT.7) I GO -7 
YZ = TS ( IGO ) 

Y 1 -TStIGO+l ) 

Y2-TS(T0C+? I 
V3"TS(TG0-»3) 

X Z - I G 0 

X = (ES-YJI*tHS-Y2»*« FS-V3WIYZ-V1J / IV2-Y2»/IY2-V3»*XZ 

1 ♦ (ES-YZ>MES-Y?)*<rS-Y3)/IYl-V2)/m-V2)/IYl-Y3)*IXZNl.l 

2 ♦CES-V71»€£S-Y1 )*C FS-V 3) / IY2-V2) / IY2-Y1 »/ *YZ-Y3J • IXZ*Z.I 

3 ♦tCS-Y?l»tE5-Yll*(FS-Y2l/IY3-Y7)/IY3-YlWIY3-Y2)*IXZ«3.l 
70 ALROT=O.S*X4l.O 

R 0 T r ] 0 . • *ALROT 
CALL IN TRP ( X t TF * F « I GO , 4 I 
C«LL IN TRP ( X . XF AYE . XU C E » ICC . 4 ) 
CPEAK:FL0ATtMAXC'!)4EF-F-XUSE*9. 

I r (IPART.r 3 . ?) CPF A K.:C PEAK-10. 

VEL - 1 0. •*( 1.329044 .01 5 05 2 * CP E AK ) 

F«CT0R:IRART4241IPSRT-1) 

TFflP-.CcC5 • V EL • VEL • F ACTOR/ROT 
CALL IN T RP C X . TZ. Zt IGC .4 ) 


ALRDT = £ 

C-f FAX - r* 

\JeL- 
Ten p -T 
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ca LL INTRF<X.TS.X c I*IBO.M) 

CALL TN T RF t X. T 3 . T 3 ! » ICO . M I 
C* LL INT.RPI X » T M . TMT . ISO. M I 
X - 3 !_R 0 T 

THECCR = 1 0.* • 1 THE T »♦ THE T# » I 5 .05 OOE-1 0. 6SE5*X *8. 6221M *X* *2-1 . 50 El 1M » 

IX. *31/10 0 0. I 

den: i .010*3 z / c toar i ipaRT i /? /(/el ♦T4ECC'? p£ N/ - s ""^p ^ 

NUTS:! 

J:13 

WRITE I 6. MO) J.Z.F.TSI.TU.TMI 

ESKE«rF 3 -T 3 I 

£KURT=EM-TMI 

COME : 0 . 

SUMC=0. 

COEF :E ■* / S OR T <E.?? 3 *ES I 
WRITE (G. 100 ) 

100 F 0 RMBT( 5 H CHAN. C X . SHUPM 8 G *6 X . d HCURM .6 X .MHCI'RT I 
DO 10 ? 1 = 1 t N C H 9 N 

X : ( T L 0 9 T 1 1 C H a N ( I I - M 9 X C H I - EF I / S Q R r IES) 

FE:£XP ( - 0. 3 * X. XI 
SKEW : X » » 3 - 3 . • X 
FKUR T rX « » M - 5 . « X * X+ 3 . 

CURM=C0EF«FE*n.'*f3*5KEW/E.*EM«FKURT/?M.I 
S UME = G UM E ♦ ( C UR M - UP M fl 0 I III * * I 
GUMC:SUMC*UPM«G(I > * UPM Q !3 I T I 
X:<FL09 T (ICH9N( TI-MSXCHI- FI /SSR TIES) 

FE:EXP(- 0 . 5 *X.X 1 
r KEW:X** 3 - 3 .*X 


GUi$_H - vtMu tj ,Uuvum£s 

^3 Cp 

C.'JX\ r ^ ^ 0 .^ 


JPL Technical Memorandum 33-519 


35 



FKUR T-X* * A -G. * X* X + 1. 

CURT=Co"F «FE • ( 1 . + T3I *SKEW/6 . »Tt|I*FKUPT/2<! . > 

WRITE IG .1 01 1ICHAN! I) .UPMAG II) .CURM ,CUR T 
101 FORMAT(I5.3F 10.0) 

1 05 C ON TINU r 

ELEG^SUME /SUMC 

WRITE 16 «1 10) ELS9.ESKEW .EKUR T £ LSQ. : A/CUV.ccWnV 

110 F ORMA TIJF10.B) - ff P 

RE TURN 
END 

3F0R.IS FMOE 

FUNCTION F MCE ( I UP . UPK A G ) 

DIMENSION V SR ( 1 7 ) .UPMAG I 1 ) 

COMMON /NOP/ NCHAN /ICHAN/ ICHAN 120) /MAXCH/ MAXCH 
00 1 T-l. NCHflN 

1 VARt I) =UPMAS !I> ‘FLOAT! IICHANII) -MAXCH) ♦ «IUP 1 
FMOE :SMP1 V AR. 1 .NCHAN) 

RETURN 

END 

3F0R.IS INTRP.INTRP 

SUBROUTINE IN TRP < X . V AP . Y . TGO .KN T ) 

C 

C INTERPOLATION 9Y BESSELS FORMULA IV ON PAGE 77 OF SCARBOROUGH. 

C IGO IS INDEX OF FIRST OF A VALUES OF VAR TO BE USED. 

C DELTA X : H = 1. 

C INPUT X. GET Y OUT. 

C 

DIMENSION VflR(l) 
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YM-tfAR(IGO) 


Y Z = VAR ( IGO+ 1 > 

DYM=YZ-YM 

IFTKNT.GT.2) 00 TO 3 
U^X-FLOA 7 (IGO) 

Y-YM+U*OYM 
RE TURN 

3 Y 1 r tf flR l IG.O* 2 ) 

OY?:Yl-YZ 
00YM-OY Z-DYM 

U- X -FLOS TtIGO+1 1 

TFJKNT.GT. 3) GO TO 4 

Y-YZ + U*(DYM + DYZ)/2 «+U*U*DDYM/2 • 

RETURN 

4 Y2-VAR(TGO*3> 

DY1-Y2-Y1 
DDYZrOYl -0YZ 
OODYM-ODYZ-OCYM 

Y = YZ + U»TYZ-*U» (U-I . > • (ODYM»ODYZ 3/4 .♦ IU - .5 ) »U» IU-1 . 3 *DDDYM/G . 
RETURN 

END 

5F0R.TS SMPiIMP 

FUNC T I ON OMP ( Y . V] . M2 ) 

DIMENSION Y t 2 3 ) 

IODD = 0 
NZ :M2 
GUM-0 » 

IF (MOD (M2-M1 »2 3 .EQ .0 ) GO TO I 
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V f M2 ♦ 1 I -0 


N 2 - M 2 + I 

1 00 1C I -M 1 • N2 
1000 = 100 D ♦ 1 

IF (I. £9. Ml .OR. 1. E9.N2J GO 70 1 
IF 1KC0 (TOOD t 2 1 .EQ.G 1 GO TC 2 
SUM:3UM*2. * Y t I ) 

GO TO 10 

2 SUMrSUM*«l. « YCI1 
GO TO 10 

'A SUM =SUM« V ( I ) 
in CONTINUE 
SMP = SUM / ?« 

RF TURN 
£ NO 
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TABLE 1. RESULTS OF FIT BY LEAST SQUARES AND HERMITE POLYNOMIALS 


# 


particle 

V 

1 

-1.6 

proton 

373.8 

2 

-2.6 

proton 

395.6 

3 

-2.5 

alpha 

391.0 

4 

-0.2 

proton 

385.6 

5 

-0.1 

alpha 

386.3 

6 

-4.7 

proton 

322.4 

7 

-5.2 

proton 

408.8 

8 

-0.7 

proton 

406.2 

9 

-4.5 

proton 

417.0 

10 

-1.3 

proton 

388.5 

11 

-2.1 

proton 

404.0 

12 

-4.7 

proton 

506.7 

13 

-6.0 

proton 

458.5 

14 

-0.8 

proton 

580.4 

15 

5.1 

proton 

617.3 

16 

r— 1 

r— < 
1 

proton 

515.0 

17 

0.4 

proton 

542.6 

18 

8.2 

proton 

544.2 

19 

2.1 

proton 

531.2 

20 

1.7 

proton 

513.0 

21 

3.0 

proton 

440.4 

22 

4.3 

proton 

589.8 

23 

5.5 

proton 

460.9 

24 

2.7 

proton 

437.0 

25 

9.2 

proton 

452.3 

26 

3.3 

proton 

439.7 

27 

1.8 

proton 

442.4 

28 

2.6 

proton 

421.2 




km/sec 


T 

n 

LS 

error 

HP 

error 

87.9 

4.1 

.0004 

.0071 

56.9 

3.1 

.0019 

.0047 

411.8 

0.4 

.0040 

.0505 

53.1 

3.4 

.0074 

.0841 

355.8 

0.2 

.0042 

.0110 

17.0 

2.8 

.0015 

.0059 

189.9 

8.9 

.0123 

.0238 

71.5 

3.5 

.0120 

.0579 

48.7 

3.6 

.0021 

.0058 

42.9 

3.4 

.0003 

.0017 

86.4 

3.4 

.0032 

.0147 

36.0 

5.7 

.0012 

.0081 

51.4 

3.2 

.0019 

.0632 

181.7 

5.6 

.0054 

.0326 

67.2 

12.5 

.0035 

.0422 

33.2 

3.3 

.0091 

.0217 

72.0 

1.7 

.0105 

.0567 

1033.9 

10.0 

.0085 

.0429 

111.2 

3.8 

.0026 

.0149 

90.2 

2.2 

.0025 

.0027 

284.9 

3.3 

.0093 

.0547 

42.9 

0.6 

.0107 

.4510 

107.9 

10.0 

.0035 

.0243 

69.4 

3.4 

.0084 

.0116 

122.0 

4.7 

.0076 

.0193 

225.2 

9.1 

.0025 

.0241 

99.2 

4.2 

.0068 

.0978 

37.6 

2.3 

.0006 

.0791 


1 ® I/ 3 

10 K cm 


eskew 

ekurt 

Av/v 

AT/T 

.01 

-.01 

.000 

.007 

.00 

-.01 

-.002 

.033 

-.12 

.20 

.003 

-.094 

-.03 

.08 

-.008 

-.250 

.08 

.16 

-.002 

.072 

.01 

-.03 

.000 

-.065 

.07 

.06 

-.001 

-.167 

.22 

-.21 

.007 

-.110 

.02 

.02 

.002 

.142 

.02 

-.03 

.001 

-.023 

-.02 

-.06 

-.006 

.150 

.03 

.03 

.000 

-.067 

.07 

.20 

-.001 

-.066 

-.16 

.19 

-.006 

-.228 

.09 

.26 

.000 

-.155 

-.09 

.19 

.002 

-.018 

-.02 

.10 

-.004 

-.224 

-.06 

-.03 

-.019 

-.391 

-.04 

.06 

-.001 

.089 

-.01 

-.08 

.001 

.004 

.06 

.06 

.001 

-.051 

.01 

-.02 

.001 

-.230 

-.12 

.10 

-.002 

-.077 

.09 

.19 

-.002 

.049 

.03 

.02 

-.000 

-.029 

.07 

-.09 

.007 

-.095 

.09 

.01 

.001 

.095 

.03 

-.04 

-.001 

-.064 


Number of points 
in spectrum 


.024 8 

.032 7 

.000 7 

-.206 5 

.000 8 

.036 5 

-.034 9 

-.029 5 

.083 5 

.029 6 

.176 7 

.070 5 

.031 5 

-.125 4 

-.104 5 

-.091 4 

-.118 5 

-.039 6 

.026 7 

.000 7 

.000 7 

-.200 3 

-.020 7 

.029 7 

-.021 7 

-.099 9 

.072 4 

-.044 5 


uo 

xO 



TABLE 2 


CHARACTERISTICS OF CURVED- PLATE ANALYZER ELECTROMETER 


Sea le 

Most Sensitive 

Middle 

Least Sensitive 

Current for 
Digital output 
= 1 , amp 

10“ 14 

0.512 

10" 12 

0.512 

10- 10 

0.512 

Current for 
Digital output 
= 511, amp 

5.12 X 10 

5_JL1_ x 10 - 9 
5.12 

5U1 , 10 -, 

Rise time 
(10 to 907„) 
seconds 

0.25 

0.025 

Rate Limited 
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Fig. 1. Logic diagram for the calculation of solar-wind parameters 














THE SUN IS NOMINALLY IN THE +x DIRECTION . 
NORMALLY INCIDENT PARTICLES HAVE 8 = 4 > = 0. 
THE xy PLANE IS THE ANALYSIS PLANE. 

Fig. 2. Definitions of 0 and <j> 


42 


JPL Technical Memorandum 33-519 


RELATIVE EFFECTIVE APERTURE = g 2 


CHANNEL CENTERS 

1-2 1-1 I 1+1 1+2 



E v 

RELATIVE ENERGY/CHARGE = — a —5 

*1 T 


Fig. 3. Relative effective aperture of the curved-plate analyzer 
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Fig. 4. The relation between x*, C^ 1 , Cj(£), and 
the channel number 
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CURRENT (ARBITRARY UNITS) 



CHANNEL NUMBER 

Fig. 6. The effect of a large EKURT 
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CURRENT, TELEMETRY UNITS 



CHANNEL 


Fig. 7. The fit of spectrum # 1 
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300 


200 


100 


9( 



:hannel 

i where |at/t| > 0. 20 


















G = GROUND 

S = SUPPRESSOR = -120 V 

M= 100 TO 11,000 V SQUARE WAVE, 2461 .5 Hz 

C = COLLECTORS 

Fig. 11. OGO 5 Faraday cup detector. (Left) Scale drawing of the OGO 5 
Faraday cups. (Right) Arrangement of the four current col- 
lectors in the OGO 5 Faraday cups and the impact location 
(shaded areas) for a beam of particles entering the instrument 
from a direction corresponding to the angles a and (3. 


JPL Technical Memorandum 33-519 


51 


A] + A 2 + Ag — 40 



Fig. 12. Values of the angles a and (3 as a function of the relative currents 
Aj, A^, and A 3 to collectors Cj, C£, and C 3 , respectively. The 
spacing of points in this figure indicate the uncertainties in a and 
P arising from the digitization of the current measurements for 
the particular case Aj + A 3 + A 3 = 40 telemetry units. 
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Fig. 13. Sample curved-plate analyzer spectra 
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